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SUMMARY 


Stresses were calculated for finite-width orthotropic laminates with a 
circular hole and remote uniaxial loading using a two-dimensional finite- 
element analysis with both uniform stress and uniform displacement boundary 
conditions. Five different laminates were analyzed: quasi-isotropic 

[ 0^/± 45^/90®] g, 0®, 90^, [0^/90^] g, and [±45^1 g. Computed results are pre- 
sented for selected combinations of hole-diameter — sheet-width ratio d/w and 
length- to- width ratio L/w. 

For small L/w. values, the stress-concentration factors were sig- 

nificantly different for the uniform stress and uniform displacement boundary 
conditions. Typically, for the uniform stress condition, the K^n values were 
much larger than for the infinite-strip reference condition; however, for the 
uniform displacement condition, they were only slightly smaller than for this 
reference. For long strips, the differences due to boundary conditions were 
small . 

The results for long strips are also presented as width-correction fac- 
tors, which relate the maximum stress in finite-width laminates to the corre- 
sponding infinite-sheet value. For d/w ^ 0.33, these width-correction factors 
are nearly equal for all five laminates and may be approximated, with a maximum 
error of less than 5 percent, by using the quasi- isotropic case. This approxi- 
mation is expected to apply to most practical laminates, provided L/w is suf- 
ficiently large. 


INTRODUCTION 

The stresses at a fastener hole are typically characterized by using an 
appropriate elastic stress-concentration factor. For metals, the necessary 
stress-concentration factors for finite-size, isotropic sheets with holes are 
usually available from the literature. However, for composites very little is 
known about stress concentrations for finite-size orthotropic laminates. The 
purpose of this paper is to explore this general problem area for composites 
and to present results for a broad range of variables. 

The study was based on a two-dimensional finite-element analysis of an 
orthotropic laminate with a circular hole and uniaxial loading. Five graphite/ 
epoxy laminates which covered a wide range of properties were selected for 
the analysis; these laminates were quasi- isotropic [ 0^/145^/90^1 g, 0^, 90^, 
[0^/90^] g, and [±45^1 g. Finite-element models were developed for a wide range 
of hole-diameter — sheet-width ratios d/w and length- to-width ratios L/w. 
These models were then cinalyzed for each laminate using both uniform stress 
and uniform displacement boundary conditions at the ends of the model. 



Results are presented in terms of stress distributions near the hole and 
as elastic stress-concentration factors for specific combinations of d/w and 
L/w for each laminate. 


SYMBOLS 


d 


Exx 


Eyy 

^xy 

^tg 

Ktn 


L 


R 


r,9 


Sn 

w 


x,y 

^xy 

c'ee 

('^ee)max 


diameter of hole, m 
Young's modulus in x-direction, MPa 
Young's modulus in y-direction, MPa 
shear modulus, MPa 

gross-section stress-concentration factor, (^90)max/Sg 

net-section stress-concentration factor, C^99)max/®n 

infinite-sheet stress-concentration factor, ('^99)max/®g 

length of sheet, m 

radius of hole, m 

polar coordinates, m and deg 

gross-section nominal stress, MPa 

net-section nominal stress, MPa 

width of sheet, m 

Cartesian coordinates, m 

Poisson's ratio 

stress component in 9-direction, MPa 
maximum oqq , MPa 


FINITE-ELEMENT ANALYSIS 

Figure 1 shows the general problem analyzed in this study. The range of 
sheet sizes and configurations was generated by using different combinations of 
d/w and L/w. Specifically, seven values of d/w were used over the range 
0.05 = d/w i 0.91, and L/w was taken as 1, 2, and 10. 

Throughout this study, the hole radius was assumed to equal unity. Accord- 
ingly, the range of d/w represented various sheet widths for a single hole 
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size, rather than a fixed width with various hole sizes. The results would be 
the same for either choice of assumptions. 

Table I presents the elastic constants for the five laminates analyzed in 
this study. As previously mentioned, these laminates were selected to cover a 
broad range of composite properties. The 0® laminate properties (ref. 1) were 
used with lamination theory to calculate the elastic constants for the other 
laminates. Although the [0^/±45^/9C^] s laminate was called a quasi-isotropic 
laminate, its calculated elastic constants satisfy all the criteria for isot- 
ropy. Consequently, the results for this laminate were compared with isotropic 
solutions from the literature. 

Figure 2 shows the finite-element models. The wide-sheet model shown in 
figure 2(a) is typical of those analyzed. The hole radius was unaltered and 
the desired values of d/w and L/w were obtained by removing elements from 
the model width and length. The narrow-strip model (fig. 2(b)) was the only 
exception to this procedure. For this case, where d/w = 0.91 , additional 
refinement was needed in the narrowest portion of the model. The corresponding 
model for L/w = 10 was generated by adding elements to the model length in 
figure 2(b). 

Two types of elements were used. For the wide-sheet model in figure 2(a), 
linear-strain (6-node) elements comprised the region near the hole (3R x 5R) ; 
uniform-strain (3-node) elements comprised the remainder of the model. The 
narrow-strip model in figure 2(b) consisted entirely of linear-strain elements. 
These models were analyzed using the BEND finite-element computer program 
described in reference 2. 

As previously mentioned, two different types of loading were considered. 
For the first case, uniform x-axis tension was specified as a boundary condi- 
tion. For the second case, uniform x-direction displacement was the specified 
boundary condition along the end of the model. The applied load corresponding 
to this uniform displacement case was calculated by summing nodal forces on the 
end of the model. For comparison, the computed local (near the hole) stresses 
for each case were normalized by the corresponding net-section nominal stress. 


DISCUSSION OF RESULTS 

Results are presented in two formats - as stress distributions near the 
hole and as stress-concentration factors. The stress distributions are briefly 
discussed to introduce the stress-concentration factors. As previously men- 
tioned, the ■’.aminates, as well as the sheet configurations, were selected to 
provide a wide range of behavior. 


Stress Distributions 

Stress distributions are presented for each laminate, using L/w = 10 and 
a range of d/w ratios. For this large L/w value, the sheet was assumed to 
approximate an infinitely long strip with uniform remote stress. 
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Figure 3 introduces the stress distributions for the [ 0°/±45°/90®] g 
quasi- isotropic laminate. Figure 3(a) shows . 099 stresses along the hole 
boundary and figure 3(b) shows this stress along the y-axis (9 = 90°). The 
099 stresses were normalized by the net-section nominal stress Sn* For sim- 
plicity, these figures show results for only three nonzero d/w values. In 
figures 3(a) and 3(b), the dashed curves represent the infinite-sheet solution 
where d/w = 0 (ref. 3). Figure 3(a) also shows the stress-concentration fac- 
tor for the infinite-sheet case. 

Figure 4 presents results for the 0^ laminate. These stress distribu- 
tions have larger maximum values (^00) max larger gradients than in the 

quasi-isotropic case. (Compare figs. 3 and 4.) These differences are probably 
caused by the large ratio of longitudinal- to-shear stiffnesses for the 0^ lami- 
nate. (See table I.) Again, the infinite-sheet results (dashed curves) were 
calculated using equations from reference 3 and the material properties from 
table I. 

The stress distributions for the 90^ laminate are plotted in figure 5. 

In contrast to the previous two cases, the largest stress magnitude is associ- 
ated with the compressive stress at 0 = 0^, rather than the tensile stress at 
6 = 90^. As expected, the [0^/90®]s results in figure 6 are approximately 
equal to the averages of those shown for the 0® and 90^ cases. 

The stress distributions for the [± 45^13 laminate have some unusual fea- 
tures. Figure 7(a) shows that (c^00)max occurs near B = 60® for d/w ^ 0.5, 
rather than at the usual location on the transverse axis (0 = 90®). Also, the 
location for this peak stress depends on the d/w ratio. Furthermore, fig- 
ure 7 (b) shows that the Oqq away from the hole may be slightly larger than 
at the hole boundary (r/R = 1.0). 


Stress-Concentration Factors 

Figure 8 shows results for the quasi- isotropic laminate; net-section 
stress-concentration factors K^n plotted against d/w for L/w = 10. 

For comparison, this figure also includes results from the literature on iso- 
tropic strips with a circular hole and uniform stress (refs. 4 and 5). Present 
results and those from the literature agree closely over the range of finite 
widths. Furthermore, the curve drawn through these results appears to agree 
with K = 3. 

oo 

Figure 9 presents all the results for the quasi-isotropic laminate. 

For each L/w value, two curves are plotted - one for the uniform stress 
boundary condition and one for the uniform displacement boundary condition. 
Results for the stress boundary condition were much more sensitive to L/w 
than were the uniform displacement results. Below L/w = 2, the stress bound- 
ary condition resulted in K^n values much higher than the uniform displace- 
ment results; at L/w ^ 2, the results were about the same. Numerical results 
are given in table II. 

The Ktn results for the 0® laminate are presented in figure 1 0 and 
table III. For L/w = 10, stress and displacement boundary conditions 
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produced the same K^n values. Therefore^ the L/w = 10 curve is a good 
approximation to the infinite-strip case, as assumed previously. The curves 
for L/w = 1 and 2 lie above or below this infinite-strip curve , depending 
on whether the sheet is loaded by the stress or displacement boundary condi- 
tion. The differences between these curves are much larger than those shown 
previously for the quasi- isotropic laminate. 

The 0^ laminate was also used to verify the finite-element procedures. 

This case was believed to be the most severe of those considered because it had 
the largest stress gradients. The results for the L/w = 10, d/w = 0.05 case 
were compared with the infinite-sheet stress-concentration factor, = 6.430, 

from reference 3. For this comparison, the computed was converted to 

gross-section stress-concentration factor K^-g. The uniform displacement and 
uniform stress conditions yielded K^g values of 6.471 and 6.517, respec- 
tively, only 0.6 percent and 1.3 percent larger than K^. 

Figure 11 and table IV show results for the 90^ laminate. For this 
extreme case, the K^n values are smaller than for the two previous laminates. 
Also, the curves for different L/w values are rather closely grouped. This 
indicates a weak L/w influence, which was attributed to the low extensional 
modulus for the 90^ laminate. 

The results for the [0^/90^]g laminate are presented in figure 12 and 
table V. The L/w influence in figure 1 2 is larger than for the isotropic 
case but, as expected, it is intermediate to the extreme behavior displayed 
by the 0^ and 90^ laminates. The K^n values for the [ 0^/90^] s laminate are 
nearly averages of those for the 0^ and 90^ laminates. 

The values for the [±45^]g laminate, presented in figure 13 and 

table VI, differ somewhat from previous results. For d/w ^ 0.50, the K^n 
results of figure 13 are based on values of (‘^00)max located near 0 = 60^ 
on the hole boundary. These K^n values are shown by solid symbols. For 
very narrow strips (d/w > 0.50), the computed K^n values are based on the 
usual location of 0 = 90^ for (O00)niax* Note that even for L/w = 1, the 
curve for the uniform displacement case is quite close to the infinite-strip 
(L/w = 10) case, in contrast to the behavior for the uniform stress case. 
Therefore, the uniform displacement case provides the better approximation to 
the infinite-sheet case. This conclusion also applies to the other laminates 
analyzed. 

Figure 1 4 is a comparison of the infinite-strip curves for the five lami- 
nates. For d/w = 0, the K^n values range from 2.48 to 6.43, compared with 
3.0 for the quasi- isotropic laminate. However, for d/w = 1, all five extrap- 
olated curves converge to a narrow range of about 1.9 to 2.3. Comparison of 
the 0® and 90^ curves shows that the d/w influence on is quite differ- 

ent for the two laminates. The 0® curve has K^n values from 6.43 to 2.3, a 
range larger than 60 percent of the maximum contrast, d/w has less 

influence for the 90^ laminate; its K^n curve shows from 2.48 to 1.9, a 

range of only about 20 percent of maximum. 

To show how d/w influences (c^60)max' results from figure 14 were 

replotted in figure 15. For comparison, the values of (c^60)niax this 
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figure were normalized by corresponding infinite-sheet values of maximum stress 
SgK^ for each laminate. Therefore, the curves in figure 15 can be interpreted 
as width corrections which relate the finite-width (c^00)max values to the 
corresponding infinite-sheet values. 

In figure 15, all five curves emanate from (^99 )max/^g^cx) " ^ diverge 

for larger values of d/w. However, these curves diverge gradually, despite 
the wide range of orthotropic laminate properties that they represent. Sur- 
prisingly, for d/w of less than 0.33, all five curves could be approximated 
by the isotropic curve, with a maximum error of less than 5 percent. This is 
believed to be significant because the present range of laminate properties and 
the range d/w 2 0.33 are expected to cover most practical cases. 


CONCLUDING REMARKS 

Stresses were calculated for finite-width orthotropic composite sheets 
with a circular hole and remote uniaxial loading using a two-dimensional 
finite-element analysis. Five different laminates were analyzed; quasi- 
isotropic [0^/±45^/90^1 s, 0^, 90^, [ 0^/90^] s, and [±45^1 s» The effects of 
sheet width and length and the differences between uniform stress and uniform 
displacement boundary conditions were studied. Stress-concentration factors 
are presented for various combinations of hole-diameter — ^sheet-width ratio d/w 
and length- to- width ratio, L/w. 

For small values of L/w, uniform stress and uniform displacement bound- 
ary conditions produced significantly different stress-concentration factors, 
especially for the 0^ laminate. Typically, results for the uniform stress 
boundary condition were much larger than for the infinite-strip reference case; 
in contrast, for uniform applied displacement the results were only slightly 
smaller than this reference. For large L/w values, differences due to bound- 
ary conditions were small. 

Anisotropy had a strong influence on the stress-concentration factors for 
long strips with small d/w values. At d/w = 0, the K^n values ranged from 
6.43 for the 0^ laminate to 2.48 for the 90^ laminate, compared with 3.0 for 
the quasi-isotropic [ 0^/±45^/90^] s laminate. By contrast, for d/w = 1, all 
extrapolated curves converged to a narrow range of about 1.9 to 2.3. 

The results for long strips were also presented as width-correction fac- 
tors, relating the maximum stress in finite-width laminates to its correspond- 
ing infinite- sheet value. For d/w ^ 0.33, these width- cor recti on factors were 
nearly equal for all five laminates and may be approximated by using the quasi- 
isotropic case, with a maximum error of less than 5 percent. This approxima- 
tion is expected to apply to most practical laminates, provided L/w is 
sufficiently large. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 18, 1979 
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TABLE I.- LAMINATE CONSTANTS 


Elastic constants 


Laminate 


Quasi- isotropic 
[0°/±45°/90°] s 

0 ° 

90° 

t0°/90°]s 

[±45°ls 




57 890 

146 900 
10 890 
79 500 
22 250 


Eyy, MPa 
57 890 

10 890 
146 900 
79 500 
22 250 


Gjjy/ MPa 
22 090 

6 412 
6 412 
6 412 
37 770 


''xy 

0.310 

0.380 

0.028 

0.052 

0.735 


TABLE II.- STRESS-CXJNCENTRATION FACTORS FOR THE QUASI- ISOTROPIC LAMINATE 


L/w 

Loading 

(a) 


s 

1 

D 

2 

S 

D 

1 n 

! S 

1 u 

D 


Diameter-to-width ratio, d/w 
0.05 0.10 0.20 0.33 0.50 0.67 0.91 

2.865 2.770 2.676 2.759 3.199 4.451 

2.850 2.714 2.468 2.219 2.008 1 .916 


2.854 


29 2.520 


2.854 I 2.728 I 2.517 




2.134 
2.1 83 2.117 

2.326 I 2.180 2.114 


6 2.326 


2.114 2.040 


®S denotes uniform stress boundary condition; D denotes uniform 
displacement boundary condition. 





































TABLE III.- STRESS-CONCENTRATION FACTORS FOR THE 0° LAMINATE 



®S denotes uniform stress boundary condition; D denotes uniform 
displacement boundary condition. 


TABLE IV.- STRESS-CONCENTRATION FACTORS FOR THE 90° LAMINATE 



®S denotes uniform stress boundary condition; D denotes uniform 
displacement boundary condition. 



TABLE V.- STRESS-CONCENTRATION FACTORS FOR THE [0O/90°]s LAMINATE 


L/w 

Loading 

(a) 

1 

s 

D 

2 

S 

D 

10 

S 

" i 


0.05 0.10 


4.626 I 4.516 


Diameter-to-width ratio 


0.33 


4.582 

4.224 ( 3.592 2.978 


. 6 ' 


6.072 

2.305 


4.592 I 4.394 J 4.061 3.738 3.367 3.014 2.400 

3.480 3.032 2.719 2.313 


4.582 4.356 3.933 3.480 


3.991 


3.165 I 2.826 2.340 

3.165 2.826 2.340 


®S denotes uniform stress boundary condition; D denotes uniform 
displacement boundary condition. 


TABLE VI.- STRESS-CONCENTRATION FACTORS FOR THE [±45°]- LAMINATE 


L/w 

Loading 

(a) 

1 

s 

D 

2 

S 

D 


Diameter -to- Width ratio 


0.05 0.10 0.20 0.33 

(b) (b) (b) (b) 


2.745 2.657 2.567 


2.741 2.639 2.480 


2.738 2.630 2.457 

2.738 2.630 2.454 


2.456 

2.456 




2.920 
1.899 

1.916 1.772 1.895 


1.904 1.767 1.895 
1.909 1.766 1.890 
1.909 1.766 1.890 


®S denotes uniform stress boundary condition; D denotes uniform 
displacement boundary condition. 

*^(‘^00)max occurred at 0 = 57°. 

*^( 000 1™=., occurred at 0 = 63°. 














































































6, deg 

(a) Stress along hole boundary. 

Figure 4.- Stresses in 0° laminate with 




1.0 1.2 1.4 1.6 1.8 2.0 

r/R 


(b) Stress along y-axis (9 = 90°) . 


iform stress boundary condition; L/w =10. 



rnrr^ 




0 30 60 90 

0, deg 


(a) Stress along hole boundary. 

Figure 6.- Stresses in (0°/90°]g laminate with 



1.0 1.2 1.4 1.6 , 1.8 2.0 

r/R . 

(b) Stress along y-axis (6 = 90°) . 
uniform stress boundary condition; L/w = 10. 



^ Isotropic sheet (Howland, ref. 4) 

^ Isotropic sheet (Koiter, ref 5) 
o Quasi-isotropic laminate (present analysis) 
□ Koo = 3.0 (ref. 3) 
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d/w 


Figure 8.- Stress-concentration factors for a quasi-isotropic laminate with 
a circular hole and uniform stress boundary cohdition; L/w = 10. 










o Uniform stress boundary condition 
/\ Uniform displacement boundary condition 
□ Kco = 4.78 (ref. 3) 



d/w 


Figure 12.- Stress-concentration factors for a [0^/90 ^] 3 laminate 

with a circular hole. 





location 


<9Q0 

0 K|-j^ from present study 
■ from ref. 3 





0 .2 .4 .6 .8 


d/w 


Figure 15.- Normalized maximum stress in orthotropic laminates with 
circular hole and uniform stress boundary condition; L/w = 10. 
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